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Abstract 

A comparative discussion of the normal form and 
action angle variable method is presented in a tuto- 
rial way. Normal forms are introduced by Lie series 
which avoid mixed variable canonical transforma- 
tions. The main interest is focused on establishing a 
third integral of motion for the transformed Hamil- 
tonian truncated at finite order of the perturbation 
parameter. In particular, for the case of the action 
angle variable scheme, the proper canonical trans- 
formations are worked out which reveal the third 
integral in consistency with the normal form. De- 
tails are discussed exemplarily for the Henon-Heiles 
Hamiltonian. The main conclusions are generalized 
to the case of n perturbed harmonic oscillators. 

1 Introduction 

The following contribution is concerned with finite 
perturbation series characterized by bounded re- 
mainders in properly chosen compact domains of 
phase space. We are, however, not interested here in 
estimating the rest terms. As has been known since 
Poincare, infinite perturbation series for Hamilto- 
nians, in general, do not converge in compact do- 
mains. Or more precisely, if convergence takes place 
according to the KAM theorem, then it generally 
occurs in an invariant subset of phase space whose 
complement is open and dense [Qj. 

Let us assume that the Hamiltonian can be 
brought into the form H = h + eV where h refers 
to n uncoupled harmonic oscillators and e is the 
smallness parameter. Then the main differences of a 
perturbative treatment by normal forms and mixed 
variable generating functions, respectively, can be 
characterized as follows. In the latter method one 



uses action angle variables / £ R n , 4> £ T n and tries 
to find a canonical transformation (/, (j)) — > (J, ^) 
which makes the transformed Hamiltonian indepen- 
dent of ip in a certain domain D of phase space. 
For an elementary introduction into this method, 
including the main ideas of the proof of the KAM 
theorem, the textbook ||| is recommended. In the 
normal form case, on the other hand, one adopts 
complex canonical variables {u u ,v v ), and tries to 
make H canonically equivalent to n harmonic oscil- 
lators given by 



h 



-jy] UJyU u Vy. 



(1) 



Both strategies fail, because terms with resonance 
denominators occur to any order, in general, which 
cannot be transformed away. Thus, even if one ac- 
cepts a finite cutoff at order e , it is not possible, 
in general, to transform a Hamiltonian into an inte- 
grable form. There seems to be at least one advan- 
tage with normal forms: they straightforwardly pro- 
vide us with a third integral. As a consequence, in 
the case of two degrees of freedom, for instance, the 
cutoff part of the Hamiltonian which is normalized 
up to order N, is integrable within the definition 
domain of the normal form transformation. In the 
case of the Henon-Heiles Hamiltonian ||] this was 
first demonstrated by @. As a further advantage, 
the normal form transformation can be carried out 
very efficiently by Lie series and thus by symbolic 
computer algebra Q. 

As a power series in e, both perturbation schemes 
should be equivalent. However, it does not seem 
to be obvious, how the third integral can be de- 
tected in the action angle variable picture. We will 
adopt an iterated transformation scheme where new 
canonical variables are introduced at each perturba- 
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tion order. Eventually, an elementary linear canon- 
ical transformation will make coming forth the ad- 
ditional integral. For demonstration, the Henon- 
Heiles Hamiltonian is considered. The results are 
generalized to the case of n perturbed harmonic os- 
cillators. 

The existence of a third integral up to order 
may be useful in Nekhoroshev-like estimates ||,||, 
||. For instance, in the three body problem an ap- 
proximate integral, in addition to energy and angu- 
lar momentum, should help to get sharper bounds 
for the remainders. We make use of this occasion to 
remark, that a former study by the present author 
H on the three body problem in celestial mechanics, 
essentially, was a failure because of explicitly and 
tacitly (eq.(77)) adopted adiabatic assumptions; the 
rigorous estimates in the Appendices C and D of 
Q, on the other hand, may be helpful elsewhere for 
similar problems. The N-body problem of celestial 
mechanics was recently examined more rigorously 
by [[h]] within the Nekhoroshev scheme. 

In the next section we will briefly introduce to the 
Lie series formulation of time dependent perturba- 
tion theory for Hamiltonian systems. Most numer- 
ical integrators are based on time series with small 
time step. Section 2 also serves to prepare the Lie 
series method for normal forms in section 3. In sec- 
tion 4 the perturbation method by action angle vari- 
ables will be discussed. In appendix A we present a 
proof of a recursive formalism for Lie series ||] , and 
in appendix B an exact canonical transformation is 
presented which brings the Henon-Heiles Hamilto- 
nian into an integrable form up to the second order 
of the perturbation parameter. 



2 Time dependent perturbation 

We assume finite dimensional systems. The Hamil- 
tonian field generates a flow which, during the time 
t, maps a phase space domain Dq of initial points 
into the domain Dt- If the field is sufficiently 
smooth in the image domains D T for r £ (0,t), 
then Dt is diffeomorph to Dq. In particular, in- 
ner points of Dq are mapped into inner points of 
Dt, and the boundaries of Dt and Dq are equally 



smooth. Moreover, by the Liouville theorem, the 
domain volumes are preserved. Chaotic behaviour, 
clearly, develops for such systems, if at all, as an 
asymptotic property. Its observability depends on 
the adopted degree of resolution. 

Lie series are tied to the time evolution of a func- 
tion f(p,q) where {p, q} = {p(t;p ,qo), p(t;qo,qo)} 
are the generalized momentum and position at time 
t of a trajectory with initial point {po,qo} £ R 2n . 
We are interested in the development of / along 
a trajectory with given initial point, and write 
therefore f(p, q) =: F(t;p ,qo) with F(0;p ,qo) = 
f(po, qo)- The Taylor expansion in the time interval 
t G (0, At) reads 



N tk 



i? W = E^ i?(fe) (i = 0) + O([Atr+ 1 ) 



fc=0 



(2) 



Defining the linear operator Lh by the Poisson 
bracket 

LHf:= ±miL_™m (3) 

^ \dp k dq k dq k dp k J 

and making use of the canonical equations, we can 
write 



\ \dp k 



df ^fdfdf. 
dt ;H \op k dq k 

= L H f. 



^Kdpk dq k dq k dp k 



(4) 



As a consequence, the Taylor expansion can be ex- 
pressed in terms of the following Lie series 



N tk 



k=0 

and in the limit ./V 
form 



(p,g)=(po,go) 



+0([At] 



N+l\ 



(5) 

oo we can write in compact 



F(t)=exp[tL H }f(p 



' c D(p,q)=(po,qo)- 



(6) 



For small enough time steps At, the time evolution 
of any dynamical variable, in particular / = p or 
/ = q, can be approximated by iterated truncated 
Lie series. A cutoff at N = 4 corresponds to a 
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fourth order Runge-Kutta integration. If H and 
/ are given analytically, the coefficients of the Lie 
series can also be determined analytically, e.g. by 
means of symbolic computer calculators. However, 
in each step symplecticity is fulfilled only up to an 
error of order [At]^" 1 " 1 . In numerical problems, it 
may be more adequate to adopt so-called symplectic 
integrators which are canonical in every step within 



the number precision of the computer, see e.g. |11]. 
The remainder of the truncated Lie series, on the 
other hand, can be rigorously expressed through (if 
/ is scalar) 



0([At] 



F N+1 (t*) 

(N + iy. 



[At] IV+i with t* G (0,At), 

(7) 



or 



with (p*,q*) = (p(t*),q(t*)). From a priori or a pos- 
teriori knowledge on the domain of (p, q) , the re- 
mainder can be estimated by upper bounds. 

3 Normal form and third integral 

For demonstration, let us consider the Henon-Heiles 
Hamiltonian 101 



H = h + eV 



with h 
V 



1 



9? 



,12 



(9) 



where {p u ,qu) £ R 2 with v = 1,2 are canonical 
variables. We assume that the variables have been 
made dimensionless, in particular H=l, which im- 
plies that the smallness parameter e is proportional 
to the square root yE of the energy of a given tra- 
jectory. With the aid of the linear canonical trans- 
formation (p u ,qu) (u u ,v u ) where 

Uu = (q» - ip u )/V2, v u = i(q u + ip u )/V2, (10) 
we obtain 



/; 
V 



-i (uivi + u 2 v 2 ); 
u 2 -iv 2 



(11) 



2^2 



(til - ivi Y 



1, 

3 (U2 



iv 2 ) 



which has the suitable form for being subject to a 
normal form transformation. 

In section 2, from the time evolution (po,Qo) 
— > (p(t),q(t)), we had derived as generator of a 
canonical transformation the operator exp[iL#]. 
Clearly, any function H{p, q) which does not ex- 
plicitly depend on time, gives rise to such a gen- 
erator. Moreover, the variable t in exp[tLjj] does 
not need to be identified as time, it can be any real 
parameter. This is seen e.g., when symplecticity 
is inferred from the Poisson bracket L p q calculated 
with respect to (po ; 9o)- It is therefore legitimate 
to adopt e as parameter of the generating function 
0. The point of view adopted here is to gener- 
ate, at any given time t, a canonical transformation 
(p(i;e = 0),<z(t;e = 0)) -> (p(t; e), q(t; e)) , which 
is parametrized with respect to the interaction pa- 
rameter e. The transformed Hamiltonian is written 
as a power series 

H(p(t;e),q(t,e))=h(p(t;0),q(t;0)) + 

oo 

^F fc (p(t;0),^;0))6 fc (12) 

k=l 

with h being the unperturbed Hamiltonian. Clearly, 
such canonical transformations can be achieved by 
means of arbitrary scalar functions. In the following 
it is convenient to write the generating function in 
the form 

X'(p(t;e),q(t;e)) : = -^x(p(i;e),g(i;e)) 

dxdp + dxdq ^ 
dp de dq de 



where we assume 



dp dh dq dh 
de dq ' de dp 



(14) 



for some scalar function h(p, q) which we do not 
need to specify. 

The canonical transformation of an arbitrary 
scalar function g is defined by the constituent equa- 
tion (we omit writing the time parameter) 



d_ 

~dl 



g(p(e),q(e))=L x ,g(p(e),q(e)), (15) 



3 



which gives rise to the power series representation 

g(p(e),q(e)) = expfeLj g(p(0), q(0)) (16) 

with X ' = x'(p(0),q(0)). 

When x, too, is expanded in a power series 

oo 

X (p(e),q(e)) = ]T e fc Xfc (p(0), <?(())), (17) 
fc=i 

then the arbitrary functions Xk will be at our dispo- 
sition to simplify the Hamiltonian coefficients H^- 
Henceforth we will write simply (p, q) for the phase 
space variables (p(e = 0),q(e = 0)). As is shown 
in Appendix A, the transformed terms, H k , can be 
determined recursively as follows ffl 



Hn 



h; L hX i + H 1 = V ; (18) 
fc-l 
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L hXk + H k = -Vk^ + Y.l^Hk-jA^) 



k = 2,3,..., where Vj is defined through the 
power series of the transformed potential, namely 
exp[eL x /]y = V ; „.,... T; <•'. 

Let us start with the term k = 1. Then we have, 
with Vq = V, 



(20) 



and we try to set H\ = with the implication 
that xi has to fulfil the relation Lh xi = V. To 
discuss, whether V is in the range of the homol- 
ogy |l2| operator Lh, we adopt the canonical vari- 
ables (10) together with the representation (1) of 
h. Furthermore, we exploit the fact that V is a lin- 
ear combination of monomials of the form U m : = 
u™ 1 D™ 2 ti™ 3 i/™ 4 with \m\ := mi+m2+m3+m4 = 3 
and rrij £ No for j = 1, 2, 3, 4. Now, each monomial 
is an eigenfunction of Lh, because 



L h U r > 



E 



dh dU r 
du v dvj, 



dh dU n 



u=l,2 

i(mi — mi + rri3 — m^)U m . (21) 



As a consequence, the set of resonance monomials 
defined by 

{U m | mi + m 3 = m 2 + m 4 = 0; mi £ Nq } (22) 



are not in the range of Lh, and therefore cannot 
be removed by the generating function xi- Clearly, 
the resonance case is possible for monomials of even 
order |m| only. Since V, according to (11), consists 
of third order terms, eq.(20) is solvable for xi with 
H\ set equal to zero. The general solution includes 
an arbitrary part of the kernel of L^ consisting of 
resonance monomials. If, as usual, this kernel part 
is set equal to zero, Xi is uniquely given by a linear 
combination of the monomials occurring in V. Fur- 
thermore, V\ = L Xl V is now determined in terms 
of 4-th order monomials. 

We examine the next iteration, which will be suf- 
ficient to reveal the general structure of the normal- 
ized Hamiltonian: 

L hX 2 + H 2 = ^V 1 + h Xl H x = ~ Vt. (23) 

Here, V\ contains both types of monomials, nonres- 
onant ones which are in the range of the operator 
Lh and resonant monomials. The latter must be 
compensated by H%, while the nonresonant terms 
are transformed away by the proper choice of X2- 
This is typical of all orders. Thus, an optimal sim- 
plification is achieved when the generating function 
is disposed of in such a way that the transformed 
Hamiltonian terms contain resonant monomials 
only. 

When this normalization is carried out up to or- 
der N, then the truncated Hamiltonian 



H (N) .._ 



N 

E 

k=0 



e k H k 



(24) 



is a constant of motion up to a rest term of the 
order e N+1 . Moreover, by the definition (22) of the 
resonance monomials and because of (21), we have 
the property 

L h HW=Q (25) 

which tells that h is in involution with and 
therefore a further constant of motion. As should be 
remarked, the remainder Rn+i in general is finite 
within properly chosen domains of phase space []l]] . 

Let E and h be the integral constants of a given 
trajectory. Going back to the original variables and 
choosing as Poincare surface of section the plane 
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q 2 = 0, one eliminates the variable p 2 from the en- 
ergy integral through 



P2 =P2{pi,qi,q2 = 0;E), 
and inserts p 2 into the third integral 
h = h{pi,q u p 2 (p 1 ,q 1 ,E)) 



(26) 



(27) 



The latter 
equation implicitly defines one-dimensional mani- 
folds M(pi, qi; E, h) = which for constant energy 
E and different values h were first plotted in refer- 
ence Q. The manifolds turned out as closed curves 
corresponding to the intersection of 2-tori with the 
Poincare plane and thus demonstrating the integra- 
bility of the approximated Henon- Heiles Hamilto- 
nian. For small enough energies E = e 2 , as is well 
known, this picture is confirmed by numerical inte- 
gration of the model. A compact symmetrized form 
of the normalized Henon-Heiles Hamiltonian can be 
found in [Ol 



4 Mixed variable generating 
function 

In terms of action angle variables (/„, fa) <— (p u ,Qu) 
defined by 



p v = \/2I v cos(fa); q v = \J2I V sin(^), (28) 

v = 1,2, the Henon-Heiles Hamiltonian reads 

H = h + eV; h = h+I 2 , 

V = V(h,I 2 ,fa,fa). (29) 

In order to reveal a third integral in the truncated 
part of the perturbatively transformed Hamilto- 
nian, we stepwise introduce generating functions as 
follows 

(Ji,j 2 , fa, fa) = Jifa + J2h + 

e n S^(J u J 2 ,fa,fa), (30) 

n = 1,2,.., which implicitly define new canoni- 
cal torus variables I 2 , fa, fa) —> {Ji,J2,fa,fa) 
through the relations 

Iv = J» + e n ^—; fa = fa + e n ^^. (31) 



In the first step, one substitutes the old action 
variables in terms of the new ones as usual to obtain 



tfM = Ji + J 2 



-irr- + "m — 



4" 



dSM T as® 

eV(J\ + e——, J 2 + e—— 



(32) 



Now we try to remove the potential term V to first 
order in e by choosing the Fourier components of 
s = Y,ni,n2ezSn 1 n 2 (Jl,J2)exp(inifa + in 2 fa) as 
follows 



S^ n2 (Ji,J 2 ) = i 



niwi + n 2 l0 2 



(33) 



Here, with the unperturbed oscillator freqencies 
uj\ = ll>2 = 1, this is possible, because resonance 
components of V with n\ + n 2 = do not exist. 
With this, the transformed Hamiltonian reads 

ff« = j 1 + j 2 + e V( 2 )(Ji, J 2 , fa,fa;e) (34) 

where, due to the elimination of the old angle vari- 
ables fa in terms of fa, the potential now is 
an infinite power series in e. 

Proceeding to second order, with the canonical 
transformation ( J\, J 2 ,fa,ip 2 ) — > ( J\, J 2 ,ipi,ip 2 ) de- 
fined by ( Ji, J 2 , fa, fa), we obtain the trans- 
formed Hamiltonian (omitting the tilde, for simplic- 

ity) 



Ji + J2 + e 2 h 2 (J 1 ,J 2 ) + 
e 2 R {2 \fa,J 2 ,fa-fa) + 

^V^{J x ,J 2 ,fa,fa;e) 



where 

h 2 {J\, J 2 ) 
/■2tt 



1 



(2vr)2 

H / * dfa V®{J x ,J 2 ,fa,fa;e 
Jo 



(35) 

(36) 
0). 



dfa 



dJ v 



It is important to realize that the angle depen- 
dence of the resonance term R^> is special and given 
through the difference fa — fa, because it contains 
only Fourier components with n\ + n 2 = 0. 
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As a consequence, if we truncate at second order 
in e, we can apply the linear canonical transforma- 
tion Ji, j 2 , i>2 — ► J{, J^'ii witn 

J{ := Ji + J 2 ; J 2 := J 2 ; ip[:=ipi; ip' 2 : = V>2 - V>i 

(37) 

to arrive at an effectively one-dimensional Hamilto- 
nian with J{ being a constant of motion. With re- 
spect to the remaining degree of freedom, (J 2 ,i/j 2 ), 
it is standard to achieve the integrable form, see e.g. 
14]. The corresponding canonical transformation is 
given in Appendix B. 

In order to see that this property continues to 
higher orders, it will be sufficient to go one pertur- 
bative step further. With the aid of the generating 
function F^ 3 \j\, J 2 , ipi, V'2) we obtain in terms of 
mixed variables 



Ji + h 



+ 



h 2 {J x J 2 ) + 0{e 3 



R {2) (Jij2,^i-^ 2 )+0(e 3 ) 

y( 3 )(J 1 ,J 2 ,Vi ) ^ 2 ;e = o) + 



0(e) 



+ 



(38) 



In the last bracket the remainder of order e stems 
from the expansion of V^(*;e) as a power series 
in e. The decisive point is that, by the chosen e- 
dependence of the generating functions (30), the 
frequencies remain unrenormalized. This is also the 
case in the normal form method. As a consequence, 
we have the same resonance condition n\ + n 2 = 0. 
Taking into account that ih y = ib v + 0(e 3 ), we ob- 
tain the third order transformed Hamiltonian (once 
more omitting the tilde) in the form 

= Ji + J 2 + e 2 h 2 (J u J 2 )+e 3 h 3 (J u J 2 ) + 
e 2 rt 2 \j 1 ,J 2 ,4> 1 -ip 2 ) + 



i R^(J 1 ,J 2 ,ij 1 -^ 2 ) + 
1 ^ (4) (Ji,J 2 ,V'i,V'2;e) 



(39) 



with the resonance terms RW,RW depending on 
the angle difference as claimed. This property, ob- 
viously, carries to the higher orders, and thus leads 



to an integrable truncated Henon-Heiles Hamilto- 
nian in agreement with the normal form. 

In every perturbation step one has to keep track 
of the definition domain of the new action variables. 
For instance, if the original variable I\ is defined 
in the positive interval [0, d%], then by (31) J\ + 

e ng S (n) JQfo is 

confined to the same domain. As a 
consequence, we have to restrict J\ to J\ S [0, d\ — 
5*} where 



5* = e n max 



ds (n) 



if 6* < di; (40) 



<£ie[0,27r] d(pi 

If 5* > d\, then the transformation is ill defined. 

5 Generalization 

The above reasoning can be immediately extended 
to the case of n perturbed harmonic oscilla- 
tors. We first discuss the normal form method. 
With the generalized multi- index notation U m := 
u™ 1 v™ 1 . . .u™ n v™ n , the eigenvalue relation (21) be- 
comes 



^ / Oh OV" Oh or 



iU m J2u„(m u -m' v ), (41) 



which gives rise to the resonance monomials 



{u r - 



m, 



0; m u £ N }. (42) 



As is remarked, even if the frequencies are all ra- 
tionally independent, there are possible resonances 
with m u = m' v for v =l,2,..n. Since, by the normal 
form method, the truncated normalized Hamilto- 
nian consists of resonance monomials only, 
eq.(41) implies the commutation of the Poisson 
bracket, namely H^O = 0, and thus establishes 
a third integral h in addition to the energy . 

In the action angle variable picture the Hamilto- 
nian (1) reads 



h = wi/i + lo 2 1 2 + ...uj n I n . 



(43) 
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It is convenient to consider the angle variables (j) v 
not on the re-torus but on the half open n-cube 
with 4> u € [0, 2?r). Adopting new canonical vari- 
ables (Ji, ... J n , ip n ) defined by the scaling (we 
assume, for simplicity, that all u> v 7^ 0) 

Iv = Jv/u) v \ 4> v = i> v LLi v 2/ = 1,2, ...re, (44) 

we obtain h = J\ + J-z + ..J n , and the Fourier rep- 
resentations 

(/ii,/j 2 ,../i„)eZ" 

exp(i H\U}\ij)\ + ... + i (inOJnipn); ip u G [0,2tt/uj u ). 

(45) 

From the resonance conditions + /i2 + ••• + /% = 0, 
we may eliminate e.g. [i\ = —[12 — ■■■ — in the 
phases of the resonance terms of the transformed 
Hamiltonian with the result that these terms de- 
pend on the following re — 1 differences only 

(46) 

Now, after the elementary canonical transformation 
(J v ,ip„) -> (Jl,ip' v ) with 

n 

(Jj,#): = (E^JuM 

(4,<): = {Ju,ih,-ifo.wi/v v ), (47) 

v = 2,3, ..re, the resonance terms do not depend 
on the new angle variable Therefore J[ is a 
constant of motion of the truncated transformed 
Hamiltonian in consistency with the normal form 
method. 

As a final remark, the Hamiltonian of the three- 
body problem in celestial mechanics (and straight- 
forwardly also the A-body case) can be expressed in 
terms of suitable action angle variables which avoid 
(chart dependent) singularities at small inclinations 
and eccentricities, see e.g. |J. It would be inter- 
esting to find out, whether a third integral can be 
worked out in a finite order perturbation procedure. 
This may be helpful in estimating upper bounds 
over finite time intervals of the order of the age of 
the planetary system. 
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Appendix A: Recursive Lie series 

We prove here the recursion relations (18) and (19) 
in a different way as compared with reference |i~5| ]. 
First we show that the coefficients of an arbitrary 
function g(p(e),q(e)) = J2k 9k £ k can be determined 
by the following recursive system Q 

9o = g(p,q); 9-n = for n = 1, 2, ..; 

n 

9n = ^~ L Xj9n-j for n = 1,2,.., (48) 

where the expansion coefficients gt and \k 
have to be taken at the point (p, q) := 
(p(e = 0),q(e = 0)). To show that (15) follows 
from the recursion system, we multiply with e n and 
sum over re 

00 00 n n 

E en 9n = 90 + E ^ E j L x 3 9n-j • (49) 

n=0 n=l j=l 

On the left hand side we have g. Differentiating 
with respect to e, transforming the double sum on 
the right hand side and making use of the fact that 
L x is linear in x, we obtain 

, 00 n 

71=1 J = l 

00 00 

= E e " _1 E-?' L X,5n-i 
n=l 3=1 

00 00 

m=0 3=1 

00 00 

= E^^XiE^m 
3=1 m=0 
00 

= y2e J ~ 1 3 L x39 = Ld 2 cg = L x ,g, (50) 

J = l 

which is (15) as was claimed. 
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For the final step we transform the Hamiltonian 
as follows 



exp[eL x ,]H= £ H k e k 

fc=0,l,.. 



(51) 



and on the other hand 



exp[e L x >] H = exp[e L x >] (h + eV) 

= E h ^ k + E v ^ k - ( 52 ) 

fc=0,l,.. fe=l,2.. 

Comparing coefficients we obtain 

Hq = h = h; H k = h k + for k = 1, 2, .. 

(53) 

Making use of the recursion formulas (48), we can 
write 

H h = E f L xi hk-j + E 7^7^ Vk-i- 3 (54) 

for k = 2, 3, ... Taking out the summand j = k from 
the first sum, and combining the remaining sums, 
we find 



fc-i 



H k - L Xk h + Y^ j,L Xj ihk-j + — -Vfc_i_j 



A; 



fc— l ■ / I \ 

= ^ + E i L x, K-j + [i + jpii^-i-i) 

= ^x fe ^+^^-i + E ^ L x, fafc-j + ^fc-i-i) (55) 

where the second term of the last equation is a con- 
sequence of the recursive system for the coefficients 
Vj. The term in the last bracket is just H k ^\. In 
view of the commutator property 

L Xk h = ~ L h Xk (56) 
we arrive at the desired recursion system (19) 

i fc-i . 

LhXk + H k = ^^ifc-i + E T^Xj ^fc-ji k = 2,3, .. 



Appendix B: Integrable second 
order form of the Henon-Heiles 
Hamiltonian 

We start from the transformed Hamiltonian (35), 
neglect the remainder V ^ , and write at first the re- 

(2) 

suiting Hamiltonian Hl runc in terms of action angle 
variables as defined in (37). We will abbreviate the 
constant of motion J{ by J. When the action angle 
variables (J^V^) are expressed by (28) in terms of 
cartesian symplectic magnitudes (p, q) , we can write 
after some efforts 



H 



(2) 



7 



7 



trunc 



J + e 2 [ J + —Jq 2 q\v + q% 

1 96 48 96 v n 

(58) 

This is one-degree of freedom Hamiltonian which 
can be brought into integrable form in a standard 
way, see e.g. [16]. The corresponding, exact canon- 
ical transformation from (p, q) to action angle vari- 
able (I, <£) is found as 



P 



-2cos($) 



(J - 1)1 



J-2sin(&)y/I(J-iy 
2T — J 

J>I. (59) 



J-28m($)y/I(J-I) 
With this we achieve the integrable form 

#£n C = </+^(l4/ 2 -14/J + J 2 

in terms of the action variables J and /. 
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